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NaaBi is a Dirac semimetal with protected nodes that may be sensitive to the breaking of time- 
reversal invariance in a magnetic field B. We report experiments which reveal that both the conduc¬ 
tivity and resistivity tensors exhibit robust anomalies in B. The resistivity p xx is B-linear up to 35 
T, while the Hall angle exhibits an unusual profile approaching a step-function. The conductivities 
cr xx and a xy share identical power-law dependences at large B. We propose that these significant 
deviations from conventional transport result from an unusual sensitivity of the transport lifetime 
to B. Comparison with CdsAs 2 is made. 


PACS numbers: 72.15.Gd,72.20.My,75.47.De,72.80.Vp 

In Dirac semimetals - the analogs of graphene in three- 
dimensional material - the bulk Dirac node is protected 
against gap-formation via hybridization. The iridate py- 
rochlores were initially predicted [I] to have protected 
nodes, but crystal growth has been problematical. Re¬ 
cently, Young et al. pj identified a class of materials in 
which time-reversal invariance (TRI) leads to node pro¬ 
tection when the nodes occur at high-symmetry points 
(the time-reversal invariant momenta or TRIM). Subse¬ 
quently, Wang et al. proposed that crystalline symmetry 
can protect Dirac nodes even when they occur away from 
TRIM. From band calculations, they identified NaaBi [2 
and CdaAs 2 [4j as Dirac senrimetals. Photoemission 0- 
i 3 and scanning tunneling microscopy ,TO] recently con¬ 
firmed that bulk Dirac nodes exist in both NaaBi and 
CdaAs 2 - Many groups HMES] predict that, in Dirac 
and Weyl semimetals, charge pumping associated with 
the chiral anomaly can be observed in an intense mag¬ 
netic field B. We report experiments on NaaBi showing 
that, even at low B , the breaking of TRI leads to ro¬ 
bust anomalies in the conductivity tensor which appears 
to originate from a strongly B-dependent transport life¬ 
time T tr (B). Among the anomalies are a robust B-linear 
magnetoresistence and an unusual step-like field profile of 
the Hall angle tan 0. The quantum oscillations show that 
Fermi surface (FS) is comprised of two distinct valleys, 
i.e. the Fermi energy lies in the non-trivial gap-inverted 
regime. 

NaaBi single crystals were crystallized from the Na- 
rich compositions (90 and 95%) tuned to preclude the 
formation of the superconductor NaBi as an impurity 
phase p3 GTS] (growth details are published elsewhere). 
The crystal structure was confirmed by X-ray diffraction 
(the lattice structure is sketched in Fig. 0c:). The deep- 
purple crystals grow with the largest facets normal to 
the c-axis (001) . To avoid deterioration of the crystals 
(which fully oxidize within 30 s of exposure to air), we 
attached contacts with Ag epoxy to the crystals inside 
an Ar glove-box and then covered them with oil before 
transferring to the cryostat (see Supplement). The resis¬ 
tivity profile p vs. T is metallic (Fig. UN with residual 
values ranging from 1.72 to 87 pUlcm. The Hall resistiv¬ 
ity Pyx is n-type and strictly B-linear (Fig. [I]), with a 


nearly T-independent Hall coefficient Rh = Pyx/B (x||I 
and z||c, where I is the current). Batch B and C sam¬ 
ples were measured without post annealing, while G1 was 
post-annealed for 1 month. Table SI in the Supplement 
lists the transport quantities measured in 8 samples. 

Figure 03 plots the MR curves in Sample B6 for se¬ 
lected 9 (the tilt angle between c and the field H = B/// 0 , 
with po the vacuum permeability). As shown by the fan 
pattern, the MR ratio p xx {B)/p xx (0) decreases rapidly 
as H is tilted into the a-b plane (9 —»• 90°). 

All the samples display prominent Shubnikov de Haas 
(SdH) oscillations in p xx (B), from which the Fermi Sur¬ 
face cross-section Sf and the Fermi wavevector kp are 
determined. In addition, we have measured the de Haas 
van Alven (dHvA) oscillations using torque magnetom- 
etry. Figure [2]A shows the oscillatory component of the 
magnetization together with a fit to the Lifshitz-Kosevich 
(LK) expression. From fits to the dHvA amplitudes vs. 
1/H and T (Panels B and C), we determine the effec¬ 
tive mass m*, the Fermi velocity vf and the quantum 
lifetime tq. The carrier density n = g v kp/3iT 2 (with the 
valley and spin degeneracies g v and g s both equal to 2) 
ranges from 2.6-4.1xl0 19 cm -3 , consistent with the Hall 
effect (Table SI). In addition, by varying 9 in the MR, 
we verify that the FS cross section Sf is nearly spherical 
(Fig. §')). 

Because of unintentional doping from vacancies, the 
Fermi energy Ep is in the conduction band, as implied 
by the n-type sign of the Hall resistivity p yx . The band 
calculations [3] predicted the existence of two Dirac nodes 
centered at (0,0, ±kp) caused by gap inversion (sketch 
in Fig. |2]D). As Ep rises in the conduction band, the 
two Dirac cones merge into a single cone when Ep ex¬ 
ceeds the Lifshitz-transition energy Ep. Recent ARPES 
experiments have confirmed the predicted dispersion and 
measured kp to be 0.095 A 1 0 and 0.10 A 1 [3] (see 
also |19j). An important issue is whether Ep in our sam¬ 
ples lies above or below Ep. As shown in Table SI, our 
measured k F ranges from 0.073 to 0.084 A . As these 
values are smaller than kp, we conclude that Ep lies be¬ 
low Ep. The FS is comprised of two distinct valleys on 
opposite sides of T, i.e. g v = 2. A persistent feature is 
a weak beating pattern in the SdH oscillations. Figure 
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& shows the SdH traces from 5 samples. The Fourier 
spectra of the oscillations reveal two frequencies f± and 
/ 2 corresponding to two values of Sp differing by ~ 16%. 
The beating suggests that the orbits may reflect quantum 
interference between the two orbits. We hope to explore 
the beating in high-mobility crystals. 

The mobility p of each sample is directly measured 
from the field profile of the Hall conductivity a xy . 
As shown in Fig. cjxy(B) has the characteristic 

dispersion-resonance profile produced by cyclotron mo¬ 
tion of the carriers. By the Bloch-Boltzmann theory, the 
extrema in a xy (B) occur at the peak fields ±S M , with 
1 /Bp = p. With increasing p. from sample BIO (p = 
21,640 cm 2 /Vs) to B6 (39,250 cm 2 /Vs) and G1 (91,000 
cm 2 /Vs), the peak field B M systematically decreases. The 
variation in p strongly influences the Hall-angle profile 
(see below). Using p = evpT tr /hkp , we find that the 
transport lifetime r* r exceeds tq by a ratio R T = 10-20 
( R t is expected to exceed 1 since tq reflects broadening 
due to all scattering processes). In Cd 3 As 2 , R T is as large 
as 10 4 [I6]. 

In the relaxation-time ansatz , the Boltzmann equation 
describing changes to the distribution function /k caused 
by an electric field E is expressed as m 


eE • v 


dfi 

dE k 


+ ev x B 


dffk 

<9k 


5k 

T~tr 
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where e is the elemental charge and cjk = /k — /k> with 
/£ the Fermi-Dirac function. and velocity v are, re¬ 
spectively, the energy and velocity at state k. The ansatz 
yields the conductivity tensor <jjj, with 


cr xx = nep/D, a xy = nep 2 B/D, (2) 

where D = 1 + (pB) 2 . From Eq. [2j the ratio <r xy /a xx = 
pB 1 which is the Hall angle tan 9, is linear in B. By con¬ 
trast, the resistivity p xx is H-independent because the 
Hall electric-field E y exactly balances the Lorentz force. 
Significantly, cr xx ~ 1 /B 2 decreases much faster than 
a xy ~ 1 /B when pB 1. These standard predictions 
assume that r tr (hence p) is a constant independent of B. 
In conventional metals and semimetals in the impurity¬ 
scattering regime (elastic scattering), this assumption is 
firmly established; the predicted trends are a cornerstone 
of semiclassical transport. 

In NaaBi, however, the observed field dependences of 
the diagonal elements a xx and p xx disagree in an essential 
way from the standard predictions (only p yx (B) appears 
conventional). As we noted in Fig. |TJ the resistivity in¬ 
creases linearly with B instead of saturating. A B-linear 
MR is rare in conventional conductors (see Abrikosov’s 
comments |20| ; we exclude metals with open orbits m- 
However, in materials with unusual topological phases, a 
B-linear MR is increasingly encountered eh To per¬ 
suade ourselves that the B-linear MR is pervasive in 
NasBi, we have investigated 8 samples (Table SI). Figure 
0A. shows that the B-linear MR is a very robust feature 
in NasBi. Across the samples, the MR ratio (measured 


at 15 T) increases from ~ 14 in B10, to 163 in G1 (the 
sample with the highest p). In B11, we show that the 
B-linear profile extends to 35 T with no evidence of de¬ 
viation. 

A second dramatic anomaly is seen in the Hall-angle 
profile. In Fig. [4j3, we compare tan 9 measured in four 
samples with increasing p, B10, B12, B6 and Gl. As 
shown, tan 9 initially rises very rapidly in weak B at a 
rate dictated by the mobility, but saturates to a plateau 
value at large B. Whereas the saturation is gradual in 
the samples with low mobility (B10 and B12), the rise be¬ 
comes abrupt in higher-mobility samples (B6 and Gl). In 
Gl, especially, the profile resembles a step-function pro¬ 
file. tan 9 assumes a virtually B-independent value from 
B = 0.5 to 15 T instead of increasing linearly with B. 
Since tan 9 = pB, the simplest interpretation of the step- 
function profile is that, starting in weak B, the transport 
lifetime varies with B as 


Ttr ~ 1 IB. (3) 

The merit of Eq. is that it also accounts for the B- 
linear MR profile, i.e. p xx = l/neT tr ~ B. Interestingly, 
with Eq. [3j the high-held B dependence of a xx is reduced 
by one power of B to cr xx (B) ~ 1/B (Eq. [2]), but leaves 
that of a xy ~ 1/B unchanged because p cancels out at 
large B. Hence both <j xx and a xy vary as 1/B at large 
B , consistent with the step-profile of tan 9. 

To verify this, we plot the B dependences of a xx and 
<j xy in log-log scale for the two high-mobility samples B6 
and Gl (Fig. [5]). In both samples, the two conductivities 
have the same power-law dependence B~@ above a rela¬ 
tively low B. Consistent with the behavior of tan 9, this 
occurs at B = 2 T and 0.3 T in B6 and Gl, respectively. 
The measured value of j3 is 1.0 in B6, but is slightly larger 
(1.15) in Gl. 

Theoretically, in the Weyl semimetal, the Weyl nodes 
which come in pairs act as sources and sinks of Berry 
curvature (Chern flux) [TJ IT^HTol . To realize a finite 
Berry curvature H(k), TRI must be broken in the Weyl 
semimetal. The 3D Dirac semimetal may be regarded as 
the limiting case when TRI is restored. In this limit, the 
Weyl nodes coincide in k space but are prevented from 
hybridizing by crystalline symmetry. Conversely [3], one 
expects that, in a Dirac semimetal, the breaking of TRI 
by an applied B will render D(k) finite, and lead to es¬ 
sential changes in its FS best detected by transport ex¬ 
periments. 

To date, theoretical attention has focused on the 
charge-pumping of the chiral current from left to right- 
moving branches [MIH [22]. However, as shown by the 
results here, the essential changes to the FS topology 
induced by B should skew a xx and p xx away from con¬ 
ventional behavior even in weak B. The results for NasBi 
(and Cd 3 As 2 [16j ) invite a systematic re-examination of 
transport in Dirac semimetals. We compare the two ma¬ 
terials in the Supplement. 

* * * 
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FIG. 1: (Color online) Magnetotransport in NaaBi. Panel A: The zero-field resistivity p and the Hall coefficient Rh vs. T 
(measured with H||c). Inset shows the crystals sealed in a vial. The largest facet is normal to c. Panel (B) shows the Hall 
resistivity p yx vs. B measured at 2 K in B6. Panel (C): The H -linear magnetoresistance in Sample B6 measured at 2 I< at 
selected tilt angles 9 to c. The MR ratio is largest at 9 = 0° (and 180°). (B = /ioH with po the vacuum permeability). The 
crystal structure of NaaBi is sketched in the inset (adapted from Ref. [5|). 
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FIG. 2: (Color online) Torque measurements of the de-Haas van Alven (dHvA) oscillations in NaoBi. The dHvA oscillations 
(solid curve in Panel A) can be fit well to the LK expression (dashed curve) with one period. From the damping versus H 
(Panel A) and the T dependence (Panels B and C), we obtain the Fermi surface section Sf, the effective mass m* = O.llmo 
(mo the free mass), velocity vf = 8.05x10 s m/s, and tq = 8.16xl0~ 14 s. The plot of the peak fields B m in and B max vs. the 
integers N (Panel D) yields Sf- 
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FIG. 3: (Color online) Panel A: Curves of A p xx showing modulation of the SdH amplitude. The beating pattern implies a 
small splitting of the fundamental SdH frequency. Panel B compares the Hall conductivity a xy (B) in 3 samples at 2 K. In the 
3 curves, the peak field ±S M yields the values /x = 21,640, 39,250 and 91,000 cnr 2 /Vs in B10, B6 and Cl, respectively. 
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FIG. 4: (Color online) Robust 77-linear magnetoresistance in NasBi (Panel A). In the 8 samples shown, p xx (B ) is measured 
with H||c at 2 K in all cases except in Bll (at 1.6 K). In Bll, the MR persists without observable deviation to 35 T. A general 
trend is that the MR increases with p, from ~ 14 (at 15 T) in B10, to ~ 160 in G1 (which has the highest p) (the MR in G1 
is plotted in | scale). In Panel B, the held profile of tan 0 = p yx /p xx is compared in 4 samples. As H increases, tan# rapidly 
saturates to an 17-independent value, which implies the anomalous relationship rt r ~ 1/H. In Gl, the change occurs at 0.5 T. 
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FIG. 5: (Color online) Log-log plots of a xx and a xy vs. B in G1 and B6. Consistent with Eq. [3] both quantities approach the 
same power law B~ 13 when B exceeds 0.3 and 2 T in G1 and B6, respectively. The measured /3 is 1.15 and 1.0 in G1 and B6, 
respectively. Curves for B6 are shifted vertically. 











